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ABSTRACT 

We investigate the correlations between optical and radio isophotal position angles 
for 14302 SDSS galaxies with r magnitudes brighter than 18 and which have been 
associated with extended FIRST radio sources. We identify two separate populations of 
galaxies using the colour, concentration and their principal components. Surprisingly 
strong statistical alignments are found: late-type galaxies are overwhelmingly biased 
towards a position angle differences of 0° and early- type galaxies to 90° . The late-type 
alignment can be easily understood in terms of the standard picture in which the 
radio emission is intimately related to areas of recent star-formation. In early-type 
galaxies the radio emission is expected to be driven by accretion on to a nuclear black 
hole. We argue that the observed correlation of the radio axis with the minor axis 
of the large-scale stellar distribution gives a fundamental insight into the structure of 
elliptical galaxies, for example, whether or not the nuclear kinematics are decoupled 
form the rest of the galaxy. Our results imply that the galaxies are oblate spheroids 
with their radio emission aligned with the minor axis. Remarkably the strength of 
the correlation of the radio major axis with the optical minor axis depends on radio 
loudness. Those objects with a low ratio of FIRST radio flux density to total stellar 
light show a strong minor axis correlation while the stronger radio sources do not. 
This may reflect different formation histories for the different objects and we suggest 
we may be seeing the different behaviour of slow rotating and fast rotating ellipticals. 
A simple analysis to estimate the effects of measurement errors indicates that the 
intrinsic degree of ant i- alignment in the roundest early-type galaxies may be as small 
as ±15° and a similar value is obtained for the degree of alignment in the late- type 
population. 

Key words: Galaxies: - radio galaxies - star-forming galaxies. 



1 INTRODUCTION 

There has been a long history of searching for alignments be- 
tween the structures of radio sources and those of their host 
galaxies and generally th e results have been incon c lusive (for 
examp l e, [M aekav (1971); Palima ka et al.l (|1 979); Va ltonenl 
|l983l ): lBirkinshaw fc Davie j (|l985l ): ISansom et ail (|l987h ). 
though in most investigations there appear to be a prepon- 
derance of objects where the radio elongation is more aligned 
with the optical minor axis than the major a xis. The clearest 
result was obtained bv lCondon et al.l (Il99ll ) who found that 
extended radio jets in 125 UGC galaxies were preferentially 
aligned with the optical minor axes of their hosts, with the 
effect being strongest for elliptical galaxies. Other, related 
investigations, have focused on the rel ative orientation of ra - 
dio axes and dust disks (for example. ISchmitt et al.l |2002)) 
and these again show a broad distribution with no obvious 
preferred alignment. Such studies were motivated by the de- 
sire to find a connection between radio emission mechanism 
and the geometry of the host galaxy. The historic results 



have been based on relatively small numbers of objects (~ 
100) but with the advent of deep radio and optical surveys 
like Fain t Images of the Ra dio Sky at Twenty centimetres 
(FIRST: iBecker et all (Il995l)) and the Sloan Digital Sky Sur- 
vey (SDSS: lYork et alJ i 20001 )). respectively, it is possible to 
construct samples with orders of magnitude larger numbers. 
This means going deeper in radio flux density to get better 
statistics and an unavoidable consequence is that the galax- 
ies identified with the radio sources are no longer dominated 
by the ellipticals, which were the targets of early studies, but 
are now a mixture of ellipticals and disk-dominated star- 
forming galaxies. Since the source of the emission is fun- 
damentally different in the two types of galaxy one might 
well expect their alignment properties also to be different. 
For example, it is known from studies of nearby star-forming 
galaxies that the rad io emission trac es the distribution of re- 
cent star-formation (|Condonl i|l992D and references therein; 
iMurphv et all l|200Sl) ) and therefore, assuming one has ob- 
servations with the appropriate resolution, one might expect 
to see a statistical alignment between the radio emission 
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and the starlight. On the other hand radio emission in the 
elliptical population is likely to be powered by the nuclear 
black hole and the alignment of the resulting jets/lobes will 
contain information about the nucleus of the galaxy. In the 
light of the existence of kinematic sub-structures, includ- 
ing decoupled cores, i n a significant fraction of early-type 
;alaxi e s (for example, Hallidav et al.l l|200ll ) ; lLoubser et al.l 
200St ): lKrainovic et alj (|2008h ). and the separation of such 
galaxies into two distin ct classes the fast and slow rotators 
|Emsellem et al.l (|2007l M. studying the alignment between 
the stellar population and the radio emission assumes added 
significance. 

In this paper we examine the statistics of the align- 
ments of radio and optical structures us ing radio data from 
the FIRST survey l|Becker et alj (119951 )1 which has a reso- 
lution 5" and opti cal data for the identifications from SDSS 
|York et al.l l|2000h ). Star formation, and its associated radio 
emission, tends to be concent r ated i n the inner few kpcof 
galaxies (for example , I Condon! (|l992l ); lMu"xlow et all (2005); 
IWilman et all (|2008h l. Thus one might expect such emis- 
sion to be just resolved in most FIRST radio maps of star- 
forming galaxies with low redshift (^ 0.2). In such galaxies 
optical elongations should be available in SDSS and thus the 
FIRST/SDSS combination is well-suited to the task. For the 
other major component of the mjy radio population, the jet- 
powered radio galaxies, FIRST data are also useful but for 
a more restricted set of objects because of the wide disper- 
sion in their linear sizes. The great advantage of using SDSS 
is the wealth of photometric and spectroscopic data that 
can be used to separate intrinsically different populations of 
galaxies and we exploit this extensively in what follows. 

Our motivations for this work are to learn more about 
the mechanisms behind the alignments and to elucidate the 
statistical properties of the low luminosity radio emitters 
which will dominate the populations of radio sources dis- 
covered in new deeper radio surveys like LOw Frequency 
ARray (LOFAR), Square Kilometre Array (SKA) pathfind- 
ers and the SKA itself. There is the possibility that the ex- 
tended radio emission from st ar-forming galaxie s might be 
used for weak lensing studies (jBlake et al.l (|2003f l) and it is 
interesting to investigate properties of the observed galaxy 
alignments when compared to the optical. 



2 SAMPLE SELECTION 

We use optical identifications of FIRST radio sources 
listed in the SDSS DR6 database (jAdelman-McCarthv et all 
ll2008t)1 whi c h wer e selected using the procedure defined in 
llvezic et all i|2002r i. There are 239993 such identifications 
which were extracted from the SDSS database using SQL. 
We have also extracted the photometric model magnitudes 
from SDSS (ugriz), the integrated flux density measured 
by FIRST (Sint), the position angle (PA, a defined between 
—90° and 90°, and measured from North through East) com- 
puted from the isophotal distribution in the r-band by SDSS 
along with the equivalent from FIRST, the major (a) and 
minor (b) axes measured by SDSS and FIRST, and R50 and 
Rgo measured from the Pe trosian int e nsity profile model fit- 
ted to the SDSS images l)Petrosianl (|l976h ). The quantity 
c = Rgo/R 5 c is known as the concentration; for c small the 
light distribution has an exponential profile and fits the pro- 



file of a disk galaxy, whereas for large values of c the profile 
is approximated by a de Vaucouleurs profile which is often 
used to model the light distribution of an elliptical galaxy. 

For our investigation we require galaxies and radio 
sources for which there are reliable measures of the PAs of 
extended emission. In order to produce a statistically robust 
sample we have imposed additional selection criteria. First, 
we have excluded all galaxies which are classified as "stel- 
lar" by the SDSS pipeline. Such galaxies are likely to be at 
high redshift and will be almost circular from the point of 
view of the SDSS beam, making computation of a meaning- 
ful PA impossible. It will also be difficult to compute the 
PA accurately for very faint galaxies, those with large axial 
ratios (b/a) and those where the major axes are significantly 
smaller than the beam. We have excluded all galaxies with 
r > 18, b/a > 0.8 in either SDSS or FIRST, and those with 
a < 2" in FIRSlfl This leaves a total of 14302 galaxies 
on which we have performed our analysis. The size of this 
sample is around 2 orders of magnitude larger than previous 
studies and it contains a range of galaxy types. 

We have decided to use the r-band PAs since they are 
computed from the images which typically have the highest 
signal-to- noise. In Fig. [1] we present histograms of the acute 
angle between the PAs measured in the next two most sensi- 
tive bands, i and g, for the sample with r < 18. It can be seen 
that there is very good agreement between them. Relaxation 
of the limit on r degrades this agreement, whereas imposing 
a more severe cut-off would produce an even tighter correla- 
tion between the PAs measured in different bands. We have 
also plotted the difference between the r and i band PAs 
as a function of b/a in Fig. [2] It shows that, as one would 
expect, the noise on measurement of the PAs is lower for 
objects with much better defined PAs, that is, those which 
have small axial ratios. 



3 RESULTS 

3.1 Whole sample 

In what follows we will present histograms of the galaxies 
in the sample binned according to the quantity Aq, de- 
fined to be the acute angle between the two position an- 
gles cvsdss and ofirst. The results for a histogram of 18 
bins are presented in Fig. [3] along with the individual dis- 
tributions for q:sdss and lyfirst- There is a clear excess 
of galaxies in which the major axes of the radio and opti- 
cal emission are aligned, Aa = 0°. Comparing the observed 
distribution with a uniform one in which bin contains the av- 
erage number per bin of the observed distribution, we obtain 
X 2 (— x 2 /degree of freedom) « 37. The number of degrees of 
freedom is the number of bins minus one since we leave the 
mean free. This can be increased to ~ 117 using 6 bins as 
illustrated in Fig. [4] We will use 6 bins as the standard in 
what follows. 

The histogram of osdss is compatible with a uniform 
distribution x 2 » 1.1. By contrast that for qfirst appears 
to have statistically significant peaks at ±30° and to a lesser 
extent at ±90°. The peaks are clearly systematics associated 



1 The FIRST beamwidth is 5" but deconvolved sizes with a < 1" 
are listed. 
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Figure 1. Comparison of the PA measured in different bands. 
The histograms are of the percentage of galaxies in each of the 
bins. The solid line represents the comparison between r and g, 
the dotted line that between g and i, and the dashed between i 
and r. In all three cases the histograms are peaked at zero and 
the differences in PA are less than 20° for ^96% of the galaxies. 



with the extra ction of PAs from the images produced the 
FIRST survey. IChang et all l|2004 ) have searched for cosmic 
shear signal due to weak gravitational lensing using the data 
from the FIRST survey. They pointed out a number of subtle 
systematic effects which had to be removed from the images 
in order to make an accurate measurement of cosmic shear. 
These effects are small at the level of an individual source, 
but since they are systematic they can be seen in global 
statistics such as those used to study cosmic shear. They 
also show up in the histogram presented in Fig. However, 
our interest is in position angle differences, Aa, and for sys- 
tematics in one set of angles to have a significant effect they 
would have to be correlated with systematics in the other 
set of angles. Since the SDSS angle distribution is nearly 
uniform it is unlikely that such systematics will have a sig- 
nificant effect on our results. Moreover, as we will see later, 
some of our most interesting results are that the observed 
position angle distributions depend strongly on the physical 
properties of the objects, and this cannot arise from mea- 
surement systematics. Nevertheless, to reassure ourselves of 
the robustness of our results we have performed some tests. 
We have selected a subset of FIRST sources to have a com- 
pletely flat distribution of position angles and we have also 
divided the original sample into different areas and looked 
for different systematics in different areas of sky. We still 
see similar biases to those in Fig. [3] in the position angle 
differences in the sample in which the FIRST position angle 
distribution is flat. In the second test we find that the FIRST 
position angle distribution is roughly the same independent 
of where we look in the sky. 



Figure 2. The percentage of galaxies against the PA difference 
measured in the r and u bands as a function of axial ratio b/a. 
The solid line has 0.6 < b/a < 0.8, dotted line 0.4 < b/a < 0.6 
and dashed line b/a < 0.4. The equivalent histogram for all values 
of the axial ratio is the dashed line in Fig. [l] 

3.2 Split on colour and concentration 

It is natural to split the sample based on photometric deter- 
mined parame ters such as colour an d concentration. It has 
been shown bv lStrateva et af] l|200ll ) that the colour defined 
by u — r can discriminate between different distributions in 
the bimodal g — r vs. u — g colour-colour diagram for an 
earlier SDSS data release. Based on morphological identifi- 
cations of nearby galaxies they suggest that u — r — 2.22 is 
the optimal colour separator and then go on to show this 
split can be linked directly to two populations, early-type 
galaxies (E, SO, Sa) which have u — r > 2.22 and late-type 
galaxies (Sb, Sc, Irr), with it — r < 2.22. We have examined 
the properties of the correlation between cvsdss and cvfirst 
as a function of u — r. Specifically we constructed histograms 
as a function of u — r with bin width A(u — r) = 0.5 and 
searched for obvious trends. We find that the peak in the 
histogram at Aa = 0° is apparent for low values of u — r 
(blue objects), whereas the histograms for higher values of 
u — r (red objects) are peaked at Aa = 90° (correlation 
between the SDSS major axis and the FIRST minor axis, 
or vice- versa) . In Fig. [5] we have split the sample into sub- 
samples with u — r < 3, 8933 galaxies, and u — r > 3, 
5369 galaxies. This appears to be close to optimal in terms 
of separating the two effects. For u — r < 3 we have that 
X 2 » 212, whereas \ 2 ~ 3.6 for u — r > 3. Clearly the 
observed correlation for blue half of the galaxy sample is 
extremely strong, but that for red half is also statistically 
signifi cant. The prec i se cut which we have used is interesting 
since IStrateva et al.l (|200ll ) showed that for u — r > 3 the 
sample would be dominated by elliptical morphologies (in- 
cluding SOs) with little contamination from spirals. We have 
also performed a split on concentration which as discussed in 
IStrateva et al.l l)200lT l can also crudely distinguish between 
different galaxy types. Using a similar approach to the case 
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Figure 3. Histograms of Act (bottom), osdss (middle) and 
<*first (top) using 18 bins for all 14302 galaxies. Included also 
are the Poisson errors for each bin. Note the peak at Aa = and 
those at ofirst ~ ±30° and ±90°. The histogram for osdss is 
consistent with a uniform distribution. 



of colour, we have identified c = 3 as being close to the opti- 
mal in separating the two effects already noted above. Fig. [5] 
shows histograms for c < 3, 8803 galaxies, and c > 3, 5499 
galaxies. The two samples have -y 2 ~ 215 and 2 respectively. 
The peak at Aa = 90° for c > 3 is less significant than the 
corresponding one when colour selection is used, suggesting 
that there is a little more cross-contamination of the sample 
when using concentration as a discriminator. 

On the basis of these two simple splits we are already 
led to the basic conclusions of this work: blue, less concen- 
trated galaxies, which broadly represent the disk-dominated, 
spiral population, have the optical and radio major axes cor- 
related, whereas at least some fraction of red, highly concen- 
trated galaxies, which are part of the elliptical population, 
have a correlation between the optical major and radio mi- 
nor axes. In the subsequent sections we will attempt to re- 
fine the selection process using other measured properties of 
these galaxies. The aim is to strengthen the correlations as 
quantified by the value of \ 2 and to obtain a physical insight 
as to their origin. 



Figure 4. Equivalent to Fig.[3]but for 6 bins. The peak at Aa = 
remains but those in »first are smeared out. The distribution 
of ojfirst is still not compatible with a uniform distribution. 



3.3 Principal component analysis 

The colour and concentration a r e kno wn to be correlated 
(see for example. IStrateva et a l. (2001)). In Fig. |S]we illus- 
trate this for our sample by presenting a scatter plot of u — r 
and c which clearly shows the degeneracy between these two 
properties. There is also a bimodality in the density of points 
in this plot, suggesting that there are indeed two populations 
within the sample. 

In order to investigate this quantitatively we have per- 
formed a principal component analysis (PCA) on the sam- 
ple, initially with two variables u — r and c. The two com- 
ponents which are generated by this procedure are 



Ci = 0.965c - 0.262(it - r) 
Gi = 0.262c + 0.965(tt - r) 



(1) 



Ci has the lowest standard deviation, 0.33, and that for 
C2 is 0.81. A line corresponding to Ci = 2.0, which is its 
mean, is presented in Fig. [6] If one considers the rotated 
coordinate system (C\, C2), the degeneracy obvious in Fig.[S] 
is along the C2 axis and therefore the value of C2 can be 
used to quantify the position along the degeneracy. Given 
that the observed bimodality appears to be compatible with 
the split along the degenerate line, we have created a set of 
histograms which split the sample in terms of the value of 
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Figure 5. Split of the sample by colour (bottom) and concentra- 
tion (top). For colour we split u — r < 3 (left) and u — r > 3 (right), 
and for concentration we use c < 3 (left) and c > 3 (right). The 
two right-hand plots are peaked at Act = 0°, whereas the two left 
hand plots are biased toward Aa = 90° . 



Figure 7. Various splits of the data determined on the basis of 
the PCA analysis and the axial ratio: (bottom left) C'2 < 3.5; 
(bottom right) C 2 > 3.5; (top left) C 2 > 3.5 and b/a < 0.6; (top 
right) C 2 > 3.5 and b/a > 0.6. 




Figure 6. Scatter plot of the colour u — r vs. c for the whole 
sample of 14302 galaxies. Included also are the lines C\ = 2.0 
(dashed) and C2 = 3.5 (solid). 



C2 and we find that C2 = 3.5 is close to optimal in terms of 
separating the two populations. 

In the bottom two panels of Fig. [7] we present results 
of splitting the sample using C 2 . There are 7573 galaxies 
with C2 < 3.5, the histogram is biased toward Aa = 0° and 
X » 260 improving the correlation from the simpler splits 
u — r < 3 and c < 3. When one creates a sub-sample with 
C2 > 3.5 there are 6729 galaxies, the bias is now toward 



Aa = 90°, and \ 2 « 5 again appearing to improve the pre- 
viously identified correlation. It seems therefore that C2 is 
a better discriminator of the populations responsible for the 
two observed correlations than either colour or concentra- 
tion by themselves. 

We have also considered further splits on the basis of 
the optical axial ratios, b/a, of the galaxies as measured by 
SDSS, in addition to C2. One reason for doing this is that 
for high axial ratios the position angle of elong ation becomes 
increa sing difficult to measure reliably. Also ICondon et al.l 
(1991) found that minor axis alignments were strong only 
for the rounder galaxies. Splitting by axial ratio appears 
to have little effect when C2 < 3.5 (not shown), but for 
C2 > 3.5 the situation is very different as shown in the top 
two panels of Fig.[7]which shows only objects with C2 > 3.5; 
that is the red highly concentrated objects which are pre- 
dominantly early-type galaxies. For b/a > 0.6 there are 5246 
galaxies and there appears to be a bias toward Aa = 90° 
with x 2 ~ 7, whereas for b/a < 0.6 there are 1483 galaxies 
with a (less significant) bias toward Aa = 0° with x 2 ~ 1.7. 
This suggests that a split based on axial ratio produces 
cleaner subsamples. It also shows that measurement errors 
on optical position angles are not strongly biasing our results 
since the stronger correlation is seen for the more circular- 
galaxies. It is interesting that splitting by axial ratio seems 
to separate two populations amongst these predominantly 
elliptical galaxies, one with Aa = 90° and the other with 
Aa = 0°. We explore this further below when we also split 
the populatio n in terms of their r adio loudness. 

Following llvezic et al.l l|2002f ) it is convenient to define 
an AB radio magnitude in terms of the integrated flux den- 
sity 
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Figure 8. Scatter plot of t — r vs. Ci. Included also are lines 
corresponding to C2 = 3.5 and t — r = —2.5. 



and the parameter t — r then quantifies whether the galaxy 
is relatively radio- loud or radio-quiet. To put the degree 
of radio loudness in context we calculate that t — r ~ 
—7 corresponds to an L* galaxy having a radio luminos- 
ity which would put it roughly at the dividing line be- 
tween radio sources hav ing FR1 and FR2 radio structures 
l|Ledlow fc Owenl (|l996l )V Though the galaxies in our sam- 
ple have a range of optical luminosities, the vast majority 
will have radio luminosities below the FR1/FR2 division. 
We note that since the beamsize of FIRST is only 5 arc- 
seconds, the measured radio flux density will be that as- 
sociated with the nuclear region and will sometimes be an 
under estimate of the total flux density. This fact should be 
borne in mind when interpreting the results presented here. 
It also should be made clear that our radio-quiet/loud divi- 
sion used later in this paper at t — r = —2.5 is different from 
the standard definition in two ways; we do not use the total 
radio flux density and, more importantly, we take the ratio 
of the radio emission to the total optical starlight and not 
the AGN lighi0. In Fig. [S] we present a scatter plot of the 
t — r versus C2 which appears to have an interesting mor- 
phology. One can clearly see that there are two populations 
separated by C2 = 3.5. The population with C2 < 3.5 has 
—2.5 < t — r < 1.0, whereas that with C2 > 3.5 extends over 
a much wider range —4.5 < t — r < 1.5. 

In Fig. [9] we have split the subsample with C2 > 3.5 
based on both b/a and t — r. For b/a > 0.6 and t — r> —2.5 
(radio-quiet; bottom left panel) there are 3441 galaxies. 
Their histogram is biased toward Aa = 90° with \ 2 ~ 11 
which is the statistically strongest result we have been able 
to find for the correlation between major axis of the optical 



2 Coincidently the numerical value for the separation between 
our radio-loud and radio-quiet populations is the same as that 
usually adopted to separ ate radio-loud and radio-quiet quasars 
llKellermann et al] |l99J) 



Figure 9. Histograms for the case C2 > 3.5 with various splits 
based on b/a and t—r. The two plots at the bottom have b/a > 0.6 
and the two the top have b/a < 0.6. The two on the left are radio- 
quiet with t — r > —2.5 and the two on the right are radio-loud 
with t - r < -2.5. 



image and the minor axis in the radio using just photomet- 
ric properties. Remarkably, for the radio louder population 
with b/a > 0.6 the distribution of the 1805 objects is statis- 
tically indistinguishable from uniform. Thus, amongst these 
galaxies with b/a > 0.6 we appear to have isolated two pop- 
ulations exhibiting different behaviour depending on their 
degree of radio-loudness. We think that the dichotomy be- 
tween radio-quiet and radio-loud galaxies for b/a > 0.6 is 
an intriguing result worthy of follow-up. We will discuss this 
further in Section 4. 

As we shall explain in section 3.4, we believe that the 
apparent trend towards Aa = 0° in the radio quieter sub- 
sample with b/a < 0.6 (top left of Fig. is due to con- 
tamination of the elliptical population by red star-forming 
galaxies. The b/a < 0.6 histogram for the radio louder ob- 
jects (top right of Fig. 0), which contains 382 galaxies, is 
biased toward Aa = 90°, albeit with a low x 2 « 1.6. We 
deem this to be compatible with uniformity. 



3.4 Spectroscopic identifications 

In the discussion so far we have relied on galaxy properties 
derived from photometric data. In addition to photometry 
many of the galaxies have SDSS spectra. These data have 
been analysed to compute the redshift of the object, but they 
have also been used by a number of authors to derive a range 
of galaxy properties such as t he stellar mass, age and star- 
formation rate. In particular, Kauffma nn et all (|2003l ) used 
SDSS spectroscopy to compute the stellar mass, m a teiiar, and 
age, isteiiar, from two stellar abs orption line ratios coupled 
to broadband photometry, while iBrinchmann et all (|2004l ) 
list star-formation rates (SFRs) and other parameters for 
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Figure 10. Histograms of the sample split on the basis of spec- 
troscopic data: (bottom left) all galaxies with spectroscopic data; 
(bottom middle) those classified on the basis of spectral lines from 
Brinchmann et al (2004); (bottom right) those with spectroscopic 
data, but were not classified on the basis of spectral lines; (top 
left) those classified as star-forming galaxies; (top middle) those 
classified as composite; (top right) those identified as AGN. 
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Figure 11. Scatter plot of c vs. u — r (bottom) and t — r vs. C2 
(top). Only objects whose spectra classify them as star- forming 
(triangles), composite (squares) or AGN (hexagons) are in the left 
hand plots, whereas those which cannot be classified on the basis 
of their spectra are on the right. Included also are lines presented 
in Figs. |5] and \S\ Note that the splits C 2 < 3.5 and C 2 > 3.5, 
appears to be virtually the same as the split into classified and 
unclassified, respectively. 



SDSS galaxies using similar technique^. We note that these 
calculations use the DR4 catalogue whereas we have used 
DR6. We have used the SDSS identifiers to match up the two 
catalogues. For the objects we have selected, all positions in 
DR 4 are consistent with thos e in DR6. 

iBrinchmann et all (|2004l ) also attempt a classification 
of the galaxies on the basis of their spectral lines, partic- 
ularly on the location of objects in the plane defined by 
ratios [Nil] /Ha and [OIHJ/H/3; that is where objects lie 
in the Baldwin-Phillips- Terlevich (BPT) diagram, follow- 
ing [BaldwineFaD l|l98ll ). The galaxies are divided into four 
groups: star-forming, active galactic nuclei (AGN), compos- 
ite (that is, have both star-forming and AGN properties) 
and unclassified. Before we discuss the results of our analy- 
sis using the SDSS spectroscopic data a word of caution is 
appropriate; the spectra were obtained with a 3 arcsec fibre 
bundle and, particularly for low redshift, face-on, galaxies 
the spectra may be missing some star-formation in disks. 
Some of these relatively rare objects could be labelled as 
unclassified (see below) because the spectroscopy only refers 
to a bulge of a disk galaxy. 

The sample of galaxies which we have been studying 
contains 6053 ga l axies which have entries in the tabl es of 
iKauffmann et all (|2003l ') and IBrinchmann et all i|2004l ). Of 
these 3752 were classified using their spectral line strengths 
and 2301 were not classified since they have either faint, or 
non-existent lines. This latter sample will also be of inter- 
est to us containing, as it does, mostly early-type galaxies, 
a sub-set of the red galaxies previously selected by their 
photometric properties. Of the spectroscopic sample con- 
taining 6053 galaxies, 1593 were classified as star-forming, 
1158 as composite and 1001 as dominated by an AGN. In 
Fig. [10] we have split the sample according to this classifi- 
cation scheme. We find that the sample of star-forming and 
composite galaxies are significantly biased toward Aa = 0° 
with x 2 ~ 114 and 97, respectively. The AGN are also 
biased toward Aa = 0°, but with a much lower signifi- 
cance (x 2 ~ 3.1). The objects which are unclassified spec- 
troscopically appear also to be biased, but this time towards 
Aa = 90° with x 2 « 3.8 

Although the numbers of objects, and hence x 2 , are 
lower for this spectroscopic sample, the percentage of ob- 
jects in each bin for star-forming and composite galaxies 
is similar to that for the photometric split with C2 < 3.5. 
Moreover, percentages for the unclassified galaxies are simi- 
lar to those for C'2 > 3.5, indicating a close link between the 
spectral classification and this photometric quantity. This 
is confirmed in Fig. 1111 where we have replotted Figs. [6] 
and [8] but now separating spectroscopically classified galax- 
ies and spectroscopically unclassified galaxies on the left 
and right, respectively. There are virtually no unclassified 
objects with C2 < 3.5 (4% contamination) and very few 
star- forming/composite galaxies with C2 > 3.5 (7% con- 
tamination). Furthermore, nearly all the classified galax- 
ies, including the AGN, are radio-quiet (that is, they have 
t — r > —2.5). This reinforces our view that C 2 > 3.5 is a 
good a photometric criterion for selecting passive elliptical 
galaxies. 



3 see http://www.mpa-garching.mpg.de/SDSS/DR4/Data/sfr 
catalogue.html for access to these data. 
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Figure 12. Histograms of the spectroscopically unclassified pop- 
ulation (which we believe are dominated by passive elliptical 
galaxies). The two plots at the bottom have b/a > 0.6 and the 
two at the top have b/a < 0.6. The two on the left are radio-quiet 
with t — r > —2.5 and the two on the right are radio-loud with 
t — r < —2.5. This plot is equivalent to Fig. [9] except that we 
have used the absence of a spectral classification, rather than the 
condition of C2 > 3.5 to pre-select the objects. 



In Fig. [T2]we focus on the unclassified population whose 
characteristic is that they have no identifiable lines in their 
spectra. We believe that this population consists predom- 
inantly of elliptical galaxies in which there is no evidence 
for current star-formation nor for the presence of an active 
nucleus. We will term these unclassified objects as passive 
elliptical galaxies. Since they appear, for the most part, to 
have C2 > 3.5 and have a bias toward Aa = 90°, we have 
attempted to improve the correlation by making splits into 
high and low axial ratios, and on the basis being radio- 
loud/quiet just as we did in Fig. [9] For b/a > 0.6 and 
t — r > —2.5, there are 1240 galaxies with a bias toward 
Aa = 90° and \ 2 « 7. The percentage of galaxies in each 
bin is similar to the case of the photometric sample with 
C2 > 3.5, although the value of \ 2 is somewhat lower due 
to the reduced number of galaxies in this sample. The other 
three sub-samples presented in Fig. [12] are b/a > 0.6 and 
t - r < -2.5, b/a < 0.6 and t - r > -2.5, b/a < 0.6 and 
t-r < -2.5 contain 724, 240 and 97 galaxies with \ 2 ~ 1.20, 
0.96 and 1.91 respectively. All three distributions appear to 
have no obvious bias toward a particular angle, and at least 
in the first two cases are compatible with uniformity. 

We briefly return to discuss Fig. [5] We see no evidence 
in Fig. [12] for the significant trend seen in the top left-hand 
panel of Fig. [9] The availability of spectroscopic information 
allows us to understand why there appears to be a significant 
non-uniformity for objects in this panel in terms of a mixture 
of star-forming galaxies and ellipticals. There are spectra for 
448 galaxies with C 2 > 3.5, t-r > -2.5 and b/a < 0.6. 228 
of these are classified as star-forming, composite or AGN 
and the histogram of their PA differences peaks at Aa = 0°, 



Figure 13. Histograms of splits of the population classified as 
AGN using their spectra, (bottom left) C'2 < 3.5; (bottom right) 
C 2 > 3.5; (top left) b/a < 0.6; (top right) b/a > 0.6. 



while there 240 which are not classified with a peak in their 
histogram at Aa = 90°. The statistical significance of the 
two peaks are much weaker than seen in the rest of the 
paper, particularly the one for the unclassified population, 
due to the reduced number of galaxies in the sub-sample. 
However, it seems clear that objects with C2 > 3.5 and low 
values of b/a contain a significant contamination from star- 
forming galaxies, composites and AGN. Visual inspection 
of the individual images confirms that these are indeed a 
population of red (w — r ~ 3) spirals. In fact 85% of the 
galaxies which are spectral unclassified have b/a > 0.6 (note 
that we have excluded galaxies with b/a < 0.8 from our 
analysis since their PA will be affected by noise). We can 
therefore refine our photometric criterion to select a clean 
sample of passive ellipticals to include only galaxies with 
C 2 > 3.5 and b/a > 0.6. 

The AGN population also deserves special mention. Not 
only is the bias toward Aa = 0° much weaker than for 
star-forming or for composite galaxies, they straddle the 
line C2 = 3.5 with almost equal numbers above and be- 
low, i.e . equal numbers of star-fo rming and elliptical galax- 
ies (see lKauffmann et al.l (|2008al ). In Fig. [T3]we present sub- 
samples of the AGN population split using C2 and b/a. For 
C2 < 3.5 there are 546 AGN which appear to be biased to- 
ward Aa = 0° with \ 2 ~ 4.2, whereas for C2 > 3.5 there are 
455 which appear to be compatible with a uniform distribu- 
tion. There is similar dichotomy when one splits according 
to b/a. There are 417 AGN with b/a < 0.6 which are bi- 
ased toward Aa = 0° with \ 2 ~ 5.2 and 584 with b/a > 0.6 
which are a compatible with a uniform distribution. It seems 
clear that dividing the sample according to the values of 
C2 or axial ratio sharpens the correlation toward Aa = 0° 
for AGN. Conditions on the axial ratio and C2 can be cou- 
pled together to create sub-samples with less contamination. 
There are 279 AGN with C 2 < 3.5 and b/a < 0.6. Their 
histogram has \ 2 ~ 5.2 and around 50% have Aa < 30°. 
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There are 295 AGN with C 2 > 3.5 and b/a > 0.6 whose his- 
togram is consistent with a uniform distribution. It is not 
clear how to interpret the uniform histograms: it could be 
that the samples are genuinely uniform, with no preferred 
orientation between the radio and optical axes or, for ex- 
ample, there could be two populations with roughly equal 
numbers, biased toward Aa = 0° and 90° 



3.5 Physical properties of the aligned galaxy 
populations 

So far we have used the spectroscopic information to classify 
(or not) the galaxies as star-forming, composite, AGN and 
uncl assified. However, we ha ve al ready alluded to the fact 
that lKauffmann et al.l (|2003f > and iBrinchmann et all (120041 ) 
deduce m s tellar) isteiiar and the SFR for a significant fraction 
of the galaxies we are studying. We therefore have an oppor- 
tunity to investigate the connections between these inferred 
properties of the galaxies and the correlations we have so far 
found in this paper. 

In Figs. [Ml and 1151 we have plotted the fraction of the 
galaxies for which spectroscopy is available as a function 
of tstdiar and tchar = m stc ii ar /SFR for the different spectral 
types. The quantity £ c har is the time which would be re- 
quired for the galaxy to form all its stars at the rate they 
are presently being produced. Not unexpectedly, we see that 
there is an increase in both quantities as one goes through 
the spectral types from star-forming, through composite and 
AGN, to the unclassified population. The star-forming and 
composite galaxies have t s t e iiai ~ 1 — 4Gyr and t c har typ- 
ically less than the age of the Universe. The AGN popu- 
lation have a very broad range of igtellar and i c har greater 
than the age of the Universe. The unclassified population 
have istellar > 5Gyr and a broad range of £ c har, typically 
~ 30 — 300Gyr. These properties are compatible with our 
interpretation of this unclassified population, as passive el- 
liptical galaxies with a dominant red stellar population and 
with the present day star-formation either being absent or 
involving a very small fraction of the total stellar mass. We 
draw attention to the fact that the AGN appear to be less 
prevalent amongst the galaxies with the smallest stellar ages, 
possibly suggesting that radio-loud activity tends to follow 
a burst of star-forming rather than being coeval with it. A 
conclusion that activity in general (not j ust radio-loud ac- 
tivity) follows a starburst was reached by ISchawinski et alj 
(2007) based on their study of SDSS active galaxies. A de- 
tailed discussion of the stati stics of radio-jets and act ivity 
in nearby galaxies is given in iKauffmann et al.l i|2008ah . 

The quantities t s t e ilar and £ c h ar can be used to the 
help elucidate some of the properties of the interesting sub- 
samples we have identified in the preceding sections. First 
let us consider the sub-sample with C2 > 3.5 and b/a > 0.6 
which was shown to have a statistically significant correla- 
tion with Aa = 90° when the galaxy is relatively radio- 
quiet and has a uniform distribution of the PA differences 
when the galaxy is radio-loud. For the most part the distri- 
butions of photometrically measured and spectroscopically 
inferred quantities are very similar for the radio- loud/quiet 
sub-samples. In particular the distributions of tsteilar and b/a 
are almost identical for the two sub-samples. The most vis- 
ibly different quantity is the distribution of t c har- Typically 
the distribution of m a tciiar is peaked at slightly larger values 
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Figure 14. The stellar age, t s tellari f° r the different spectral 
types: (bottom left) star-forming; (bottom right) composite; (top 
left) AGN; (top right) unclassified. Typically t s t e Uat ~ 1 — 4Gyr 
for star-forming and composite galaxies, whereas it is much larger 
(^stellar > 5Gyr) for the unclassified (passive elliptical) popula- 
tion. Interestingly the value of t s tellar spans a wide range of values 
for the AGN population. 
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Figure 15. The characteristic age, t c hau for the different spec- 
tral types: (bottom left) star-forming; (bottom right) composite; 
(top left) AGN; (top right) unclassified. Note that the scale is 
logarithmic. The star-forming and composite populations have 
*char typically less than the age of the Universe, whereas for the 
AGN and unclassified populations it is greater than the age of 
the Universe. 
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Figure 16. t cria r f° r galaxies with C2 > 3.5 and b/a > 0.6, that 
is, those in the bottom two panels of Fig. 9. The solid line is the 
radio-quiet population with t — r > —2.5 and the dashed line is 
the radio-loud population with t — r < —2.5. 
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for the radio-loud sample, whereas the SFR is peaked at a 
lower value. This combines to yield a difference in the distri- 
bution of tchar between the two samples which is presented in 
Fig. 1 16 1 The median value of t c har for the radio-quiet galax- 
ies is around a factor of 2 lower than that for the radio-loud 
galaxies. 

We have shown in Fig.[l3]that the AGN population can 
be split using the parameter C2 into a sub-sample which has 
a significant correlation toward Aa = 0° (C2 < 3.5) and 
one with an almost uniform distribution of PA differences 
(C2 > 3.5). The properties of the stellar populations of the 
AGN hosts when split on the basis of C2 are presented in 
Fig. 1171 We see that the distribution of SFR is very similar 
for the two sub-samples, but that objects with C2 > 3.5 typi- 
cally have a larger m stc n ar implying that the t c har is typically 
larger. The sample with C2 < 3.5 has t c har sharply peaked 
around the age of the Universe similar to that for the com- 
posite population, whereas the distribution of t c har is much 
broader and peaked at around 30 Gyr for C2 > 3.5 more 
in keeping with that of the unclassified population. The dis- 
tribution of isteiiar is also interesting. For C2 < 3.5 we find 
that tsteiiar is typically < 5Gyr whereas for C2 > 3.5 it is 
> 5Gyr. 

The above results suggest that there could be a corre- 
lation between tsteiiar and the parameter C2 for AGN. In 
Fig. \TE\ we plot thenar against C2. There is a clear corre- 
lation; there is reasonably small scatter about a line given 
by 

tstell a 



Gyr 



> 0.5(8C* 2 - 19). 



(3) 



There is no hint of bimodality in this distribution. It is in- 
teresting to note that star-forming and composite galaxies 
would typically reside in the bottom left of this diagram 
(tsteiiar < 5Gyr and C2 < 3.5) whereas the unclassified pop- 
ulation would reside in the top right (tsteiiar > 5Gyr and 
C2 > 3.5). The position angle differences for the AGN with 
young stellar populations appear to share the same prop- 
erties as the other galaxies in that quadrant. On the other 
hand, the older AGN appear not to share the same prop- 
erties as the other galaxies in their quadrant, though the 



Figure 17. Properties of the AGN population split on the basis of 
C2: (bottom left) M stc iiar; (bottom right) tsteiiar! (top left) SFR 
in M0yr -1 ; (top right) t c har- The solid lines are for C2 > 3.5, 
whereas C2 < 3.5 for the dashed lines. 
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Figure 18. Scatter plot of t s tellar v s. C2 for galaxies which were 
identified as being dominated by an AGN. Included also is the 
line C2 = 3.5 which has been used to split these galaxies into 
a sub-sample biased toward Aa = 0° and one where the PA 
differences are uniformly distributed. There appears to be a strong 
correlation between the two quantities. 



numbers are small and the corresponding statistical uncer- 
tainties are large. 
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3.6 Intrinsic alignment between the radio and 
optical structures 

An important question is the strength of the intrinsic align- 
ment between the radio and optical emission both for the 
star-forming galaxies and the passive ellipticals. In both 
cases the distributions of measured position angle differences 
are broadened by measurement noise. The observed variance 
is due the intrinsic effect {(Aq) 2 )i and measurement noise 
((Aa) }m, which we will assume are independent and can, 
therefore, be added in quadrature 

{(Aq) 2 ) = {(Aq) 2 ) : + <(Aa) 2 } M . (4) 

In what follows we will attempt estimate the value of 
((Aq) 2 )!' 2 based on a few simple assumptions. Our inten- 
tion is to give an idea of the level of correlation rather than 
presenting a robust statistical analysis. 

In addition to measurement noise, there can be projec- 
tion effects. These are likely to be negligible in the case of 
star-forming galaxies which can be assumed to be disk-like 
with roughly co-spatial radio and optical emission, b ut could 
be im portant in ellipticals which can be triaxial iBinnevI 
(1978). The possible effects of triaxiality on the observed 
radio /optical positio n angle distribution were simulated by 
ISansom et al.l (| 19871 ). Their simulations show that, for oblate 
spheroidal galaxies, those with jets along the minor axis 
should show the kind of alignments we observe, whereas 
other configurations lead to much less well defined position 
angle difference distributions; any triaxiality only serves to 
blur to the kind of biases we obs erve. An important recent 
contribution has been made by iPadilla fc Straussl (|2008l ) 
who have analysed the shapes of 282203 elliptical galaxies 
from SDSS. They reach the conclusion that the intrinsic 
shapes of ellipticals are consistent with them being oblate 
spheroids (with the degree of flattening increasing with de- 
creasing luminosity). The implication of ellipticals being 
oblate is that projection effects should not be significant 
and, therefore, we will ignore them 

The measurement error, quantified by {(Ao) 2 }m, can 
be reduced by focusing on galaxies which are bright both in 
the radio and in the optical and, in addition, have low values 
b/a. However, some caution is required since it is possible 
that the correlation we are interested in is intrinsically con- 
nected to the value of b/a (or even the optical or radio flux); 
for example, galaxies with low b/a might have an intrin- 
sically smaller dispersion in position angle differences. We 
have found no evidence of this in the case where the bias 
is toward Aa = 0° with all the results for C2 < 3.5 and 
star-forming galaxies being compatible with measurement 
noise being dominant in the histograms when considered as 
a function of b/a. The situation is more complicated in the 
case of C2 > 3.5 and passive ellipticals, where the strength 
of the correlation appears to improve as one increases the 
lower limit on b/a. This suggests that the true correlation 
is stronger than we have observed. We will assume that we 
can measure ((A«) 2 )m as a function of b/a for star-forming 
galaxies and that this holds for elliptical galaxies. 

Again for objects with C2 < 3.5, we have made cuts 
to the sample so as to give the strongest possible results 
in terms of the percentage of objects in each bin. We note 
that this is a somewhat different philosophy to that used 
in the rest of the paper where we have attempted to max- 



imize the x 2 - The results of this process are presented in 
Fig. [I9j We have imposed stringent cuts of r < 15, in or- 
der to improve the accuracy of determination of the optical 
PA, and Si n t > 5mJy, to do the same for the radio PA. We 
will make the assumption that for b/a < 0.5 and C2 < 3.5 
the observed correlation is purely from the intrinsic scatter 
and that measurement errors are negligible. By performing 
simulations with different Gaussian dispersions, and ignor- 
ing the very few outliers, we deduce that the dispersion seen 
in the case where b/a < 0.5 is compatible with a standard 
deviation of ((Aq) 2 )|' 2 ~ 15°. This is a conservative ap- 
proach since by not explicitly taking into account any radio 
position angles measurement errors, any intrinsic spread we 
deduce will be an over estimate of the true value. Although 
the existence of a correlation between the radio and optical 
PAs in these systems may not be unexpected, the fact that 
it is so strong is possibly a little more surprising. 

Again for objects with C2 < 3.5 we now consider the 
correlations when b/a > 0.6 (which is presented Fig. I19J1 
and b/a > 0.7. We find that the histogram with b/a > 0.6 is 
compatible with a total dispersion of ((Aa) 2 ) 1 ' 2 ~ 30°, and 
for b/a > 0.7 it is compatible with ~ 40°. Assuming that 
these galaxies are from the same parent population as those 
with b/a < 0.5 - that is they have the same intrinsic scatter 

- and that the intrinsic effect is uncorrelated with the effect 
of the axial ratio - that is they can be added in quadrature 

- we can deduce ((Aa) 2 )^ 2 . We find that it to be « 25° for 
b/a > 0.6 and ^ 37° for b/a > 0.7. 

For the C2 > 3.5 population, essentially the passive 
ellipticals, we have already shown that the correlation to- 
ward Aa = 90° is strengthened as b/a increases. The his- 
togram with b/a > 0.6 is compatible with a total dispersion 
of ~ 45° . If we assume that the measurement error arising 
from large axial ratios is the same as in the case of C2 < 3.5, 
then one can deduce that the combination of the intrinsic ef- 
fect (((Aa) 2 ), 1/2 ), and any residual effect due to projection, 
is ~ 37°. For b/a > 0.7, the histogram is compatible with a 
total dispersion of ~ 40° and therefore one can deduce that 
the intrinsic plus residual projection effects (assumed small 
if they are oblate spheroids) give a dispersion of ~ 15°. This 
is a surprisingly small value, especially given that previous 
inve stigations, albeit with much smaller numbers (for exam- 
ple, [Sansom~etjil] (|l987)), failed to produce a statistically 
significant result. Physically the likely implication of such 
a small dispersion is that we are dealing with well ordered 
systems in which the overall shape of the host galaxy dic- 
tates the nuclear accretion axis and thus the spin axis of 
the central black hole. The fact that the intrinsic dispersion 
appears to be smaller in the more circular systems is again 
unexpected. 



4 DISCUSSION 

In the previous section we have identified and refined two 
correlations between the optical and radio PAs; one is be- 
tween the two major axes (Aa = 0°) and the other is 
between the major and minor axes (Aa = 90°). In the 
simplest terms star-forming, predominantly spiral, galaxies 
have Aa = 0° and passive ellipticals have a bias towards 
Aa = 90°. We have tried to refine the division between ob- 
jects displaying the different types of behaviour and in ta- 
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Figure 19. Histograms of the percentage of objects in each bin 
when the sample is cut so as to maximize the number of objects 
with a particular correlation. The bottom two plots are selected 
to have C 2 < 3.5, r < 15 and S; nt > 5mjy. In the plot on the 
left, we have b/a < 0.5, whereas on the right we have b/a > 0.6. 
The top two plots have C 2 > 3.5, r < 15 and S- ln t > 5mjy. This 
time we plot those with b/a > 0.6 on the left and b/a > 0.7 on 
the right. 



bles[T]and[2]we present a summary of the various cuts made 
on the data based on photometric and spectroscopic prop- 
erties respectively. We find that a new composite parameter 
C 2 = 0.262c + 0.965(ii — r) can be used to discriminate bet- 
ter between different populations, with most objects with 
C 2 < 3.5 having identifiable emission lines in their spectra, 
and those with C 2 > 3.5 typically being passive ellipticals 
with no identifiable lines. Use of the C 2 parameter improves 
the statistical significance of the position angle difference 
separation. The propensity to select passive ellipticals can 
be increased by additionally requiring that b/a > 0.6. Inter- 
estingly, galaxies identified as being dominated by an AGN 
appear to straddle the line C2 = 3.5. In these, as in the 
galaxy population as a whole, AGN with younger stellar 
populations and those with low axial ratios have Aa = 0°, 
whereas those with older stars which are also more circular 
have a uniform distribution of PA differences. 

There is an obvious interpretation for the very strong 
bias toward Aa = 0° in the galaxy population as a whole 
which is that it arises from the subset of all galaxies in which 
there is significant star-formation and where the radio emis- 
sion comes, either directly or indirectly, from supernova rem- 
nants associated with the evolution of recently formed mas- 
sive stars. There is little direct evidence from the SDSS and 
FIRST data alone that the emission at the two frequencies 
is accurately co-spatial but one should note that the resolu- 
tion of the radio images is a factor of ~4 worse than that of 
the optical images. Also single colour optical images are a 
relatively inefficient way of distinguishing regions of active 
star formation from the general stellar population. Hence 
there is clearly scope for 1 arcsec resolution radio mapping 
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Table 1. Summary of the splits made on the data on the basis of 
the spectroscopic properties. N is the number of objects satisfying 
the criteria and \ 2 is th e \ 2 P er degree of freedom when compared 
to a uniform distribution. 



Split 


N 


x 2 


Bias 


All 


6053 


56 


0° 


Classified 


3752 


114 


0° 


Unclassified 


2301 


3.8 


90° 


Star-forming 


1593 


97 


0° 


Composite 


1158 


39 


0° 


AGN 


1001 


3.1 


0° 



Table 2. Summary of the splits of the data based on the classifica- 
tion using emission lines. Corresponding histograms are presented 
in Fig. 1101 We used the same notation as in tabic HI 



and more detailed optical/infrared imaging of selected SDSS 
objects to investigate in d etail the correspon dence between 
radio and optical emission. iPatel et alj (|2009h are investigat- 
ing in detail the shapes of high redshift star-forming galaxies 
and their associated emission using Hubble Deep Field data 
together with high resolution VL A/MERLIN observations. 

The bias towards Aa = 90° in the red, objects with 
high concentrations is less strong but still highly statistically 
significant. Since in these predominantly elliptical galaxies 
the emission is almost certainly powered by mass accretion 
on to a central disk/black hole system, the radio elongation 
we measure will be related to the overall spi n axis of that 
system . Our results confirm those reported bv lCondon et al.1 
l|l99ll ) who saw a strong trend for the jet axes in UGC ellipti- 
cal galaxies to be aligned with the opt ical minor axes. Sim- 
ulations done by ISansom et al.l (Il987f ) for different galaxy 
geometries show that one only gets the kind of strong mi- 
nor axis alignments that we observe for galaxies which are 
oblate spheroids. We therefore interpret our results as imply- 
ing that there is a bias for the spin axis of the central engine 
to be aligned with the minor axes of galaxies and also that 
these galaxies are predominantly oblate spheroids. That el- 
liptical galaxi es are oblate sphe r oids i s consistent with the 
conclusions o f lPadilla fc Stra uss (2008) based on an analysis 
of galaxy shapes. The radio minor axis alignment is some- 



thing one might expect to see if the galaxies were rotation- 
ally supported and the black hole accretion disk axes were 
aligned with the overall galaxy rotation axes. One should 
contrast what we see here in elliptical galaxies with what 
seems to be a much more complex situation that appears to 
be common in Seyfert galaxies where there is only a weak 
relationship between j et axis and the ax i s of h o st galaxy disk 
1 Kinney et all l|2000l ); iGallimore et al l l|200d ); iRaban et all 
( 20091 )1. 

The most remarkable thing that we find about the 
Aq = 90° bias in the elliptical galaxies is the fact that strong 
alignment only seems to be present in the radio quieter sub- 
set of these objects. We emphasise that we have not tried to 
optimize our dividing line in radio loudness to maximize the 
discrimination between the position angle properties. The 
division at t — r = —2.5 separates the samples into sub- 
samples with a ratio in numbers of around 2:1 and was not 
chosen a priori to have any particular physical significance 
for the elliptical galaxy population. However, our results sug- 
gest the division does have some physical significance and 
is consistent with there being two types of object within el- 
liptical population exhibiting quite different behaviour. We 
reiterate that the two radio populations we talk about here 
have nothing to do with the traditional FR1 /FR2 radio mor- 
phological division between that occurs at around t — r = — 7 
not t — r = —2.5. 

Based mainly on optical morphological and kinematic 
data, it has been convincingly argued that there are in- 
deed two sub-populations of elliptical galaxies; there are 
those that are fairly round, are slowly rotating and are 
generally massive, and there are those with more elliptical 
isop hotes, are fast rotatin g and are generally less massive 
fsee lEmsellem et al.l l|2007t )). These distinctions also appear 
to be relate d to w h ether or not the g alaxy has a core (see 
iFaber et all (|l997f ); iKormendv et al.l (|2009l ') and references 
therein for a comprehensive discussion of elliptical galaxy 
properties). What we have found by including radio posi- 
tion angle information suggests that such a split may be 
evident down to the scale of the central black hole. In radio- 
loud objects alignment of the central engine appears to be 
independent of the overall shape of the host galaxy, while 
in the radio-quiet population the engine knows about the 
galaxy shape. The degree of radio-l oudness is known to be 
a strong function of galaxy mass (|Best et all <|2005l )) and 
the more massive galaxies tend to be of the non-rotationally 
supported type. Thus it is very tempting to suggest a sim- 
ple picture in which the radio-quieter objects are predom- 
inantly rotationally supported with matter on all scales 
sharing a common rotation axis, while the louder ones are 
not rotationally supported and have the complex kinemat- 
ics often present i n the rounder and most massive galaxies 
|Emsellem et all i|2007l )1. However, the fact that amongst 
the radio quieter population we see a stronger minor axis 
alignment in the rounder galaxies does not have a natural 
explanation. Putting this caveat aside, in the above picture 
having complex kinematics, perhaps even a decoupled core, 
results in a more efficient radio-jet producing engine, but one 
that can point in any direction with respect to large-scale 
shape of the host galaxy. 

To investigate the radio-quiet/radio- loud dichotomy 
further we have looked for other differences in measured and 
inferred properties of the two groups using the data from 
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iKauffmann et alj l|2003t ) and lBrinchmann et all (|2004) . The 
only significant differences we can find are that median val- 
ues of the stellar mass and t c har are a factor of ~2 lower in 
the radio-quiet population. These two quantities are not in- 
dependent since the star-formation rates in the two types of 
galaxies are approximately the same and i c har is the ratio of 
the stellar mass to t 8 t e iiar- Since it is known that the fraction 
of galaxies having radio luminosities above some given limit 
increases rapidly with increasing host galaxy stellar mass 
|Best et al.l l|2005l )). the smaller average stellar mass for the 
radio-quiet galaxies is only to be expected. We have looked 
explicitly to see if the degree of alignment we see depends on 
optical luminosity (or stellar mass) and find that there is no 
obvious dependence. It is noteworthy that the lower charac- 
teristic age for the stars amongst the radio-quiet population 
argues for a more abundant supply of gas from stellar mass 
loss in these and this only serves to emphasise the point 
that the efficiency with which the hot gas available from the 
whole galaxy is converted into radio emission must be signif- 
icantly higher in the ra dio-loud objects. Of cours e cold gas 
can be accreted too (see Kauffmann ct al. (2008a)) but, be- 
ing passive elliptical galaxies, we ex pect objects in our sam- 
ple to be in the regime described by IKauffmann fc HeckmarJ 
( 2008) in which accretion is limited by the supply of hot gas 
from evolved stars. Thus the physical origin of apparent di- 
chotomy remains a mystery but is probably related in some 
way to the formation history of the galaxies since the effect 
we observe depends on both the large-scale and small-scale 
structure. In turn this could be related to the environment 
in which the galaxy has evolved. For example, the radio-loud 
objects could live in dense environments around groups and 
clusters, whereas the radio-quiet objects could live in more 
quiescent environments. This issue is presently under inves- 
tigation (Battye et al., in preparation). An alternative hy- 
pothesis could be that the intrinsic shape of the radio-loud 
galaxies is not oblate and hence projection effects hide any 
intrinsic correlation which might exist. 

The behaviour of the AGN population is also interest- 
ing. Galaxies hosting active nuclei span the full range of 
galaxy types. Active galaxies with young stellar populations 
have the radio and optical axes aligned, albeit somewhat 
more weakly than in the case of star- forming galaxies. This 
suggests that in these the radio emission is dominated by 
that from supernova remnants which are correlated with 
starlight. The situation is less clear for galaxies with an older 
stellar population hosting an AGN. These AGN appear to 
share nearly all of the properties of the spectral unclassified 
population indicating that they are in all other respects pas- 
sive ellipticals. However, there is no statistically significant 
trend for their radio emission to be aligned with the host 
galaxy minor axis (Fig. I13p . the trend we see in the ellip- 
tical population as a whole. Perhaps this is because of the 
relatively small numbers involved. 

On the basis of the apparent correlation between t s t e iiar 
and C2 for AGN and the absence of any kind of bimodality 
in Fig. [18] it might be tempting to suggest that there could 
be an evolutionary track from the the blue-extended region 
to the red-compact region, as opposed to two separate popu- 
lations, and that AGN possibly represent disk galaxies in the 
process of turning into passi ve ellipticals. Such a con clusion 
is supported by the work of ISchawinski et al.l (|2007t ). How- 
ever, we think that this is probably an over-simplification 
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since the time-scale for the activity is likely to be much 
shorter than that required to affect the transition from a 
disk-dominated system into an elliptical. 

There is one simple but quite significant consequence 
following from the fact that we identify two different align- 
ment behaviours in the general elliptical galaxy population. 
It argues for the radio emission that we observe being qui- 
escent over long time-scales. It can be seen from Fig. [11] 
that a range of t — r of 5 (a factor 100 in flux density ratio) 
encompasses virtually all the objects in the undivided popu- 
lation. Thus it would require variability of no more than an 
order of magnitude to blur the distinction between the align- 
ment properties of the radio-quiet and radio-loud objects out 
of existence. Plausible timescales on which aligned objects 
might be transformed into a misaligned one are difficult to 
assess but probably likely to be greater than the dynamical 
timescale, perhaps a few x 10 s yr. Stability of the radio lu- 
minosity on this kind of timescale would argue for a stable 
source of fuelling, most likely mass loss from the stellar pop- 
ulation consistent with the con clusion recently presented by 
iKauffmann fc Heckmanl (|2008l ). 

Finally, an original motivation for this study was to in- 
vestigate what might be found when very deep radio surveys 
are made with, for example, the SKA in order to detect 
weak lensing. Such studies are of interest because optical 
measur ements may be susc eptible to a range of systematic 
effects. IChang et all l|2004 ) have detected cosmic shear us- 
ing the FIRST survey and claim that, in principle, the ra- 
dio should be less effected by systematics. More conserva- 
tively, whatever they may be, the systematics are likely to 
be different and hence uncorrelated. Therefore, correlating 
the signal between the two wavebands will provide a cleaner 
result. Our results suggest that the radio direction of elon- 
gation of galaxies found in the very deep surveys required 
for weak lensing studies will be the same as the optical with 
a dispersion of « 15°. 



5 CONCLUSIONS 

The results we have presented show that consideration of the 
relationship between the optical and radio position angles 
can lead to a number of interesting inferences. The power of 
these results has been significantly enhanced by the rich and 
varied astronomical information available from the SDSS. 
We believe that our analysis introduces a new slant on the 
very significant knowledge of galaxy evolution which has al- 
ready come from the SDSS. We have shown that galaxies 
with young stellar populations which are forming stars at 
a significant rate predictably have aligned optical and radio 
axes with a dispersion of around 15°. Much fainter versions 
of such star-forming galaxies will be useful for future weak 
lensing studies and should complement optical studies of the 
same fields. 

Perhaps more surprisingly we have found a correlation 
between the radio major axis and the optical minor axis in 
radio-quiet passive ellipticals which, once we take into ac- 
count the range of noise and projection effects related to the 
orientation, have a similarly low dispersion to that seen in 
the star-forming gala xies. Such a co rrelation has been long 
sought (for example, iMackavl (119711 ) ) because of its poten- 
tial to tell us about the intrinsic shapes of the galaxies and 



the relationship between distribution of stars and the engine 
that generates the radio jets. Our results argue strongly in 
favour of there being a population of oblate spheroidal galax- 
ies with radio jets aligned with the stellar minor axes. We 
find, in fact, evidence for two populations of passive ellip- 
ticals. Those which are relatively radio-quiet show strong 
radio alignment with the optical minor axis while the ra- 
dio louder population has no apparent correlation between 
the optical and radio emission. We tentatively identify our 
two populations with those independently found from op- 
tical photometric and kinematic studies of ellipticals, the 
rounder non-rotating and massive galaxies, and the more 
disky, rotating and less massive ones. The rotationally sup- 
ported disky galaxies will naturally be oblate and have a 
single well-defined axis. On the other hand, the slowly ro- 
tating galaxies may have complex kinematics with the jet 
axes misaligned with the main stellar distribution. Jets may 
be bent as they propagate and triaxiality will give rise to 
projection effects. All this can go some way towards explain- 
ing the lack of radio/optical alignments in the radio louder 
population. There is clearly great scope for further detailed 
optical and radio observations of these relatively radio- weak 
elliptical galaxies. 
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